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850. 
ON LINEAR DIFFERENTIAL EQUATIONS. 


[From the Quarterly Journal of Pure and Applied Mathematics, vol. xxt. (1886), 
pp. 321—331.] 


1. THE researches of Fuchs, Thomé, Frobenius, Tannery, Floquet, and others, relate 
to linear differential equations of the form 
a™ d™— 
hiath ocd t+ +Pmy =0, 
or say E 
(Po; Pis cers Pm dz , 1)" y ca 0, 


where Po, Pis -+> Pm are rational and integral functions of the independent variable s; 
any common factor of all the functions could of course be thrown out, and it is 
therefore assumed that the functions have no common factor. It is to be throughout 
understood that # and y denote complex magnitudes, which may be regarded as points 
in the infinite plane; viz. æ, =&+%im, is the point the coordinates whereof are £, 7: 
and similarly for y. 


2. Suppose #a—a is not a factor of p, the noint #=a is in this case said to 
be an ordinary point in regard to the differential equation; and let Yə; Jı, Yo, +++, Ym-ı 


be arbitrary constants denoting the values of y, f SE. eee oe for the point «=a. 


We can from the differential equation, and the equations derived therefrom by successive 
differentiations in regard to æ, obtain the values for «=a of the subsequent differential 


2 dy dmy 2? 
coefficients dam’ qa’ e> Say these are Ym, Ymi,--»; Viz. the value of each of 
these quantities will be determined as a linear function of Yə, y, -> Ym; and we 


thus have a development of y in positive integer powers of «—a, viz. this is 


1 
Y=Ptpn@e-at7_y@—ayt..., 
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which will in fact be a sum of m determinate series multiplied by Yọ Y Yo: <e, Ym— 
respectively; say the form is 


Y =YoXot Nh + tee + Ymm, 
where Xo, X,, ..., Xm are each of them a series of positive integer powers of x—a. 
Each of these series will be convergent for values of œ sufficiently near to a, or say 
for points x within the domain of the point a; and since Yo, Yı, ..., Ym— are arbitrary, 
each series separately will satisfy the differential equation; and we have thus m 
particular integrals of the differential equation. 


3. Suppose next that w—a is a factor of p, the point «=a is in this case 
said to be a singular point in regard to the differential equation. The foregoing pro- 
dm dy agmty 
aun Zs ...3 and we 
have to consider the developments of y which belong to the neighbourhood of the 
point «=a, or say to the domain of the singular point «=a. This has to be done 
separately in regard to each of the singular points, and among these we have, it may 
be, the point œ. To decide whether œ is or is not a singular point, we may in the 
differential equation write «=1/t; and then transforming to this new variable, and 
throwing out any common factor, the coefficients of the transformed equation will be 
rational and integral functions of ¢, without any common factor; say these are 
Py, Py, ..., Pm; if t is not a factor of P,, then œ will be an ordinary point of the 
original equation; but if ¢ is a factor of P,, then œ will be a singular point of the 
original equation. 


cess of development fails, as leading to infinite values of 


4. In considering the singular point w=a, we may, it is clear, transform the 
equation to this point as origin, and it is convenient to do this; supposing it done, 
we have an equation wherein #=0 is a singular point, viz. p, contains æ as a factor, 
and we have to consider the developments of y which belong to the domain of this 
singular point s =0. 

It is convenient to change the form of the differential equation by dividing the 
whole equation by the first coefficient p, and then expanding each of the quotients 


fe D n, 2m in a series of ascending powers of æ. The new form is 


mn pi 
od +p, a A . + pny = 0, 
or say 
P(y)=(, pis sis PETA 1)" y = 0, 
where Pı, Pa «+», Pm now denote each of them a series (finite or infinite) of integer 


powers of æ, but containing only a finite number of negative powers of a, 


5. Such an equation frequently admits of a “regular integral” of the form 
y=«°H(«), where E(x) is a series of positive integer powers E, + Es + L,a?+..., 
(E, not =0, for this would imply a different value of p)*. To determine whether 


* The expression regular integral is afterwards used in a more general sense, see post, No. 10. 


50—2 
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this is so, we substitute in the differential equation for y the value in question 
æE (æ), thus obtaining a series 


NA + ON, aP-*) +... 


{where @ is a determinate positive integer depending on the negative powers of æ in 
the equation); the coefficients 0,, 2,,... are functions of p of an order not exceeding 
m, and contain also the coefficients Æ,, #,, #,,... linearly; in particular, Q, contains 
E, as a factor, say its value is = Æ,Il,. ‘The series should vanish identically. Supposing 
that Il, contains p, then we have II,=0, an equation of an order not exceeding m 
for the determination of p. For any root p=p, of this equation, Æ, remains arbitrary 
and may be taken =1; the equations 0,=0, 0,=0,... then serve to determine the 
ratios to Æ, of the remaining coefficients Æ,, Æ, ...; and we thus have the solution 
y= a (1+ Ex + Ea + ...), where p, and the coefficients have determinate values, 


6. I stop to notice a curious form of illusory solution; the assumed form of 
solution is 
y=æ(... + Esa Enst E, + Ex +...) 


the series being a double series extending both ways to infinity, or say a back-and- 
forward series; we have here a series of equations 


..O.,~=0, 0.,=0, Q,=0; Q, =0,...; 


which leave p undetermined, but determine the ratios of the several coefficients to one 
of these coefficients, say Æp; or taking this =1, we have a solution 


y =a (...4+ Hw? + Ew 4+14+ Et Ea + ...) 


where the coefficients are determinate functions of the arbitrary symbol p. Such a 
series is in general divergent for all values of the variable, and thus is altogether 


without meaning. As a simple instance, take the differential equation a -—y=0, 


which is satisfied by 
y=|. (p= 1) p= par + ap t pars ty ety 


see my paper, Cayley, Note on Riemann’s paper, “ Versuch einer allgemeinen Auffassung 
der Integration und Differentiation,” Werke, pp. 331—344; Math. Ann. t. Xvi. (1880), 
pp. 81, 82), [751]. 


7. A more general form of integral is Thomé’s “normal elementary integral,” 


y= ew (x), where w is = a finite series “+. +2 of negative powers of æ (a a 


positive integer, =2 at least). To discover ale fd a form exists, observe that, 
p , E’ (a) 


writing for a moment wy a y’, and so for the other symbols, we have A =w ia ee TE (a)? 
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where the last term 7 may be expanded as a series of positive integer powers 


of x, so that, writing iene we have 


(a — 1) eo 
acces ae 


c 
z= n 4 E +a +Q; 


we can, by introducing into the differential equation the new variable z(= A in place 


of y, obtain for z a differential equation (not a linear one) of the order m —1, and 
this should be satisfied by the series in question, viz. a series containing a finite 
—(a— 1) ca 


number of negative powers of æ; we endeavour to determine the first term = 


’ 


or say =% where the exponent a is a positive integer which is to be determined, and 
A is not =0; this is done by a well-known process; we write in the equation 
z =: and (if possible to do so) determine a (a positive integer =2 at least) by 
the condition that two or more of the terms shall have the same negative index — p 
preceding in order all the other indices, viz. the absolute value of p must be greater 
than the absolute value of any other negative index; we then equate to zero the 
whole coefficient of the term in question 2’, and thus obtain a value not =0 of A. 


And having thus obtained the first term >, we can, by assuming the form 


with indeterminate coefficients, and substituting in the equation, find the remaining 
coefficients B,..., K, L; we then have 


sie Ca G ? uk f 
giving the value ae At sae of w; instead of determining the subsequent coefficients 


from the z-equation, we may, from the original equation, writing therein y=e”Y, obtain 
an equation for Y(=e-”’y) which will be linear of the order m, and will admit of a 


regular integral Y=a?°H (a). Or we may, going on a step further with the z-equation, 
find therefrom the coefficient L which is =p. 


8. As an example, take the equation 


oy S Faye Bus 
ata aatar? 


x 
which has the elementary normal integral 


1 
z 


=e zi 
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To investigate this, writing A =z, we have Lorre, and thence the z-equation 


Substituting herein the value 2-4, we obtain the function 


AG: Be A AR 
Tant gaT gate Te 


and here the value a=2 gives the two terms s and af each with a negative 


index — 4, the absolute value whereof is greater than the absolute values of the other 
indices; we then have A*—A=0, giving a value A=1, which is not =0; and we 


have thus a first term tae E the complete value is, as it happens, the finite series 


a’ 


, ; 1 
£= -5; and z being =f , we thence obtain the integral in question, ee +e 


ii 
a m 

9. Returning to the question of the regular integrals, we may consider the 
function l 


P (a), i (lp ht: 1)” 2, 


say this is the so-called “determinirende Function,” but I will call it the Indicial 
function; this is a function of (x, p); we have therein terms arising from 


Ee a 


viz. these are equal to [p]”"a?-™, [p] "a-m, ... respectively, but these will be multiplied 
by powers of æ contained in the coefficients p,, po, ..., Pm respectively; the function 
is thus of the degree m as regards p, but the coefficient of any power of æ is of 
the degree m, or of some inferior degree, according to the terms [p]”, [p]”~,... which 
enter into the coefficient. The coefficient of the lowest power of æ in the indicial 
function may be termed the indicial coefficient, or simply the indicial; and the 
equation obtained by equating this coefficient to zero the indicial equation. By what 
precedes, the indicial is a function of p, which is of the degree m at most, but which 
may be of any inferior degree, or even be an absolute constant. Hence, considering 
a regular integral aH (æ), the index p is determined by the indicial equation, and 
to each root of this equation there corresponds in E a regular integral we (æ); 
the number of regular integrals is thus at most = 


10. Suppose, first, that the indicial equation is of the degree m, and that the 
roots of this equation are all unequal; there are, in this case, m values, say pi, Pa, ---, Pm 
of p, and each of these gives rise to a regular integral #°# (æ), viz. there will be corre- 
sponding to each value of p a series of positive integer powers Æ (æ), which will be 
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a convergent series for sufficiently small values of æ. Some of these series may be 
finite series, viz. this will be the case if the indicial equation has roots differing from 
each other by integer values (but a finite series may of course be regarded as con- 
vergent); and in the case of two or more equal roots, it is necessary to extend the 
notion of a regular integral by including under it series multiplied into positive 
_integer powers of log æ; but, with this modification, the conclusion holds good; if the 
indicial equation be of the degree m, the number of regular integrals will always be 
=m, viz. there will be m integrals involving series Æ (æ), or, it may be, E(x) (log æ}, 
where each series Æ (æ) is a convergent series for sufficiently small values of 2. 


11. But suppose the indicial equation is of an order m—vy inferior to m, and to 
fix the ideas suppose that the roots are all unequal. There are in this case m-— y 
values, SAY pı, Ps, +++» Pm—y Of p, and each of these will apparently give rise to a 
regular integral #°H(#); but it may be that in any such case the series F(a) is a 
series divergent for all values, however small, of the variable æ, and that we thus 
have in appearance only, but not in reality, a regular integral; or say the integral is 
illusory. The conclusion is that the indicial equation being of a degree m -— y inferior 
to m, the number of regular integrals is at most =m —y; but it may have any less 
value than m—y, or even there may be no regular integral. It is hardly necessary 
to remark that, if the indicial be an absolute constant, or say if the degree of the 
indicial equation be =0, then there is no value of p, and consequently no regular 


integral. 
12. By way of illustration, consider first the differential equation 


for which the degree of the indicial equation is equal to the order of the equation 
(=2); viz. the indicial function is 
=p(p—1)ar*+ par *+ a, 
and the indicial equation is thus p(p—1)+p=0, that is, p*=0, viz. there are here 
two roots, each = 0. 
It is convenient not in the first instance to write p=0, but to substitute in the 
equation the value 


y = Aw + Ba + Cart +... 
We thus find 
wr, a, WP ss. 
O=p.p—1.A, p+1.p.B, p+2.p+1.C 
+ pA, +p+1.B, +p+2.0 
+ A, +B, 
viz. we have the equations 
j pA =Q; 
(p+1¥PB+4=0, 
(p+2~C0+B=0, 
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where observe that in each equation the function of p, which multiplies the posterior 
coefficient, is of a degree superior to that which multiplies the other coefficient. The 
first equation gives p=0; A arbitrary, but say A=1; the values of B, C,... are then 


1 1 1 
——~,, Ha 9 ea eros 
(p+1) < (+IP A+ (e+ 1 (p +2P(p + 3) 
giving rise to a series which in this form, and when p is put equal to its value, 
p=0, is at once seen to be convergent (even for indefinitely large values of æ, but 
this is an accident; it would have been enough if the series had been convergent 
for sufficiently small values), viz. the series is 


£ a? a 
y=1-7,+ 7 pie get 


There is another integral 
bell ies Sheet l 2a — 3a + Ah a — 52h, at 
y= -rtr og @ + 2a — $a? + hh a — ytt., 
involving log # and a second series which is also convergent. 


13. But consider next the before-mentioned equation 


dy (8 1\dy,1.. 
a * (G-z) dated 
for which the degree of the indicial equation is equal to order —1(=1); in fact, here 
the indicial function is 

= p (p — 1) 2°? + 3px — par + a, 
and the indicial is thus =—p, giving the index p=0. There is thus at most one 
regular integral H(«), but this is in fact an illusory one. Substituting in the equation 


the value 
y= Aæ + Burr + Cort? +..., 


we find 
Ort, a, a, aP. 
0= p(p —1) 4, (p+1)pB, (p+ 2)(p +1) 6, 
— pA, —(p+1)B, -— (9+2)C, -— (9 +8)D, 
+ 3p, +3(p+1)B, 4+3(p4+2)C, 
+ A, S By Iit C, 
and we thus have the equations 
pa = 0, 


(p+ 1 A-(pt+ 1) B=0, 
(p +2" B-(p +2) C =0, 
(p +8) C —(p +8) D =0, 


where observe that in each equation the function of p, which multiplies the posterior 
coefficient, is of a degree inferior to that which multiplies the other coefficient. The 
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first equation then gives p=0; A arbitrary, or say =1: the other coefficients B, C,... 
are then found to be 


p+1, (p+1)(p+2), (p+1)(p +2) (e+3),... 


or substituting for p its value =0, the series is 
1+1e+1.2.0+1.2.32*+..., 


which is divergent for any value whatever, however small, of 2; or the integral is as. 
mentioned an illusory one. 


14. The last-mentioned differential equation is reducible, viz. we have 


elated tL IAN id eda yh 
Pa)= {aa (o a) ante} 9 (deta) (ta a)» y = Oe 
where of course the order of the factors Q, D is material. 


It is a general theorem that the indicial, belonging to the product P, is of a 
degree equal to the sum of the degrees of the indicials, belonging to the factors 
Q, D respectively; and, in fact, the indicial of P is (as was seen) of the degree 1; 
and the indicials of Q, D are of the degrees 1, 0 respectively. But, moreover, if the 
indicial of P is of a degree m—y, less than m the order of the equation (in the 
present case m=2, m—y=1), then, in order that the equation may have m —y regular 
integrals, it is a necessary and sufficient condition that P shall be decomposable into 
a product QD, where the factors Q, D (being functions of the same form as P) are 
of the orders y and m —y respectively, and where, moreover, the second factor D shall 
have an indicial of the degree m—vy, and consequently the first factor Q en indicial 
of the degree zero, In the present case, it is the second factor which has an indicial 
of degree zero; and thus the condition is not satisfied. And accordingly the equation 
ought not to have (m—y=) 1 regular integral; and we have seen that there is in 
fact no regular integral. The investigation serves to show in what sense it is that 
there is no regular integral, or generally in what sense it is that the number of 
regular integrals may be less than m—y. 


15. The theorem may, it appears to me, be stated in the more general form: the 
necessary and sufficient condition that the differential equation P(y)=0 of the order 
m, but having an indicial of the degree m—vy, shall have m—y—6 regular integrals 
is that the function P shall be a product of the form P=QMD, where the orders 
of Q, M, D are y+5—90, 0, m—y-—6 respectively, and the degrees of their indicials 
5, 0, m—y—6 respectively ; or, to denote this in a compendious manner, say 


m y+ô-0 0 m=~y-8 


my ô 0 m—y—5 


? 


@ is an integer which is not =0, and which is at most = y, for otherwise the function 
c. XII. 51 
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Q of the order y+8 -8@ could not have an indicial of the degree & If there is no 
regular integral, then m=y+6, D is a mere constant or say unity, and the formula is 


where @ is not =0, and is at most =m — y. 


16. For differential equations of the form P (y)=0 above considered, Floquet gave 
a theorem, which as afterwards remarked by him is not generally true, viz. this was 


P=ABC..., a product of linear factors of the form La Cia’ (i an integer). The 


question of decomposition is considered in my next following paper “On Linear Differ- 
ential Equations (the Theory of Decomposition),” [851], and by means of the formule 
there given (and by the process indicated ante No. 7), it could be decided whether any 
particular differential equation admits of a decomposition P=ABOC..., where the linear 


factors are functions not of the form just referred to, but of the form $ 4 30,0, in 


which the series contains only a finite number of negative powers, say this is a 
decomposition into linear regular factors. 
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